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Preface 
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processing.  Messrs.  Leo  Koestlei  ai«l  Roger  Meier,  WES,  and  Prof.  Glenn 
Rix,  Georgia  Instimte  of  Technology,  assisted  in  field  measurements  at  the 
WES  test  site.  Ms.  Jennifer  J.  Davis-Smith,  a  co-op  student  from  Mississippi 
State  University,  helped  to  prqiare  rqiort  figures.  Dr.  Mary  Ellen  Hynes  was 
die  Chief,  EESB,  and  Dr.  Albert  J.  Bush  HI  was  the  Chief,  CDAB,  during 
die  study. 

Overall  direction  at  WES  was  provided  by  Dr.  A.  G.  Franklin,  Chief, 
EEGD,  Dr.  George  M.  Hammitt  n.  Chief,  PSD,  and  Dr.  William  F. 
Marcuson  HI,  Director,  GL. 
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1  introduction 


Federal  and  state  highway  officials  are  finding  through  experience  that 
existing  design  codes  (e.g..  Applied  Technology  Council  1981  and  1983; 
American  Association  of  State  Highway  and  Transportation  Officials  1991) 
may  be  inadequate  to  mitigate  damage  produced  to  bridges  by  moderate  to 
large  earthquakes  especially  at  soft  soil  sites  (Miranda  1993).  A  number  of 
different  research  ai^  development  efforts  are  underway  at  various  institutions 
to  examine  this  problem.  One  such  research  program  ^ministered  by  the 
Federal  Highway  Administration  (FHWA),  is  aimed  at  geotechnical  aspects  of 
seismic  design  for  bridge  abutments,  embankments,  and  foundations. 

Paramount  to  calculating  the  seismic  behavior  of  a  bridge  are  the  earth* 
quake  motions  at  the  foundations.  Seismic  forces  are  ^plied  at  the  abutments 
and  at  piers  along  the  length  of  the  deck.  The  soil  conditions  at  each  founda¬ 
tion  are  likely  to  be  significantly  different,  however,  possibly  varying  from 
rock  to  soft  soil  to  rock  and  almost  always  varying  in  thickness  of  soil  below 
the  abutment  or  thickness  of  soil  below  a  pier.  Wave  passage  effects  and 
asynchronous  motion  may  also  impact  the  selection  of  appropriate  seismic 
motions.  Combinations  of  these  factors  may  produce  a  large  variation  in 
seismic  inputs  over  the  length  of  the  bridge. 

Traditional  solutions  involving  equally-distributed  seismic  forces  iq>pear  to 
be  in  question  and  more  comprehensive  dynamic  response  analyses  of  the 
complete  bridge-foundation  system  or  individual  foundation  systems  are  being 
considered.  The  more  comprehensive  analyses  allow  for  variations  in  founda¬ 
tion  and  soil  conditions  and  earthquake  wave  directivity.  Hence,  a  thorough 
understanding  of  how  to  adequately  cbaracto’ize  and  idealize  soil  stiffiiess 
(elastic  moduli)  in  the  vicinity  of  the  bridge  foundations  and  within  embank¬ 
ments  is  warranted. 

Extensive  research  and  development  of  seismic  methods  to  determine  elas¬ 
tic  moduli  has  occurred  over  the  past  decade.  These  methods  include  the 
downhole,  uphole,  crosshole,  and  the  Spectral-Analysis-of-Surface-Waves 
(SASW)  meAods.  Data  collection  and  processing  for  nearly  all  of  these 
methods  were  developed  under  the  assumption  that  the  ground  is  level  and  the 
sub-surface  layers  are  horizontal  and  extend  to  infinity.  These  assumptions 
generally  do  not  iq)ply  to  bridge  abutments  which  are  typically  placed  at  the 
crest  of  a  sloping  (iq)proadi)  embankment  fill  or  a  cut  slope. 
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The  overall  research  program,  entitled  ”In-Situ  Determination  of  Dynamic 
Soil  Properties  for  Hi^way  Structures,"  is  intraded  to  carefolly  examine  the 
effects  of  complex  geosystems  (soil-rock  systems)  and  anisotropic  states  of 
stress  on  existing  seismic  methods.  The  efforts  will  include  numerical  analy¬ 
sis,  laboratory  testing,  and  field  measurements  using  existing  techniques  and 
possible  modifications. 

The  seismic  methods  to  be  considered  as  part  of  the  program  are:  down¬ 
hole,  crosshole,  SASW,  and  seismic  tomography.  This  report  (Rq)ort  1  of  a 
series)  specifically  summarizes  the  analysis  of  the  vertical  component  of  sur¬ 
face  waves  in  sloping  ground  using  numerical  analysis  and  field  measure¬ 
ments.  Once  each  method  has  been  re-evaluated,  a  procedure  for  applying 
various  seismic  geophysical  methods  to  the  characterization  and  idealization  of 
complex  soil-rock  systems  will  be  developed. 

The  SASW  method  offers  a  unique  and  economical  alternative  means  for 
the  determination  of  stiffness  for  near-surface  conditions.  One  simplifying 
assumption  made  for  inversion  schemes  of  the  SASW  method  is  that  the 
ground  is  level  which  not  only  precludes  the  existence  of  discontinuities,  but 
also  implies  that  the  state  of  stress  is  isotropic.  The  presence  of  a  slope  defies 
this  assumption.  Therefore,  surface  waves  propagating  up,  down,  and  along  a 
slope  were  evaluated  using  the  results  of  numerical  approximations  and  field 
measurements  to  assess  the  effects  of  sloping  ground. 

To  begin,  the  SASW  method  is  described  in  Chapter  2  of  this  report.  A 
brief  description  of  the  numerical  analysis  procedures  and  validation  of  the 
computer  c^e  used  are  presented  in  Cbiqiter  3.  The  results  of  the  numerical 
analysis  which  considers  several  different  slope  angles  and  set-back  distances 
are  described  in  Chapter  4.  Field  measurements  conducted  at  a  test  site  at 
WES  are  described  and  compared  with  numerical  results  in  Chapter  S. 
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2  The  SASW  Method 


The  analysis  of  surface  waves  has  been  used  for  several  decades  to  obtain 
the  variation  of  soil  stiffness  with  d^th.  The  steady-state  surface  wave 
technique  was  developed  u'  the  late  19S0*s  (Heukelom  1958;  Jones  1958; 
Heukelom  and  Foster  1960)  for  the  analysis  of  pavement  systems.  This  tech¬ 
nology  was  ad2q)ted  by  WK  to  be  an  exploration  tool  at  soil  sites  (Fry  1963 
and  1965;  Ballard  19^).  This  technique  involves  interactively  dttermining 
foe  wavelength  at  several  frequencies  of  excitation  by  measuring  foe  distances 
from  foe  source  to  successive  standing  wave  (displacement)  peaks.  The  veloc¬ 
ity  calculated  knowing  foe  wavelength,  X,  and  frequency  of  excitation,  / : 

V*  =/X  (2.1) 


was  assumed  to  correspond  to  foe  Rayleigh  wave  velocity,  Vn,  at  an  average 
dq)fo  between  X/3  and  X/2. 

The  SASW  method  was  developed  nearly  two  decades  later  (Heisey  1982; 
Heisey,  Stokoe,  and  Meyer  1982;  Nazarian  1984;  and  Nazarian  and  Stokoe 
1985a,  1985b)  under  foe  direaion  of  Professor  Stokoe  at  foe  University  of 
Texas  at  Austin  (UT).  Since  that  time,  research,  development,  and  applica¬ 
tion  activities  have  continued  at  UT  under  the  tutelage  of  Professors  Stokoe 
and  Roesset  (e.g.,  Calderon  1985;  Nogueira  1986;  Sanchez-Salinero  1987; 
Sheu  1987;  1988;  Kang  1990;  Sykora  1993),  at  foe  University  of  Michi¬ 

gan  (e.g.,  Hiltunen  1988;  and  Gucunski  1991),  and  at  WES  (Alexander  1992) 


The  SASW  method  was  established  to  obtain  a  more  accurate  variation  of 
dastic  moduli  with  dq;>fo.  The  procedure  involves  measuring  seismic  signals 
at  foe  ground  surface  (thus  it  is  non-intrusive)  produced  by  either  a  harmonic 
or  impulsive  surface  excitation.  In  foe  original  inversion  schemes,  foe  site 
was  assumed  to  be  ideally  layered  foorizontal  and  extending  to  infinity)  and 
only  plane  waves  were  assumed  to  propagate.  Newer  inversion  schemes 
developed  at  UT  also  provide  for  foe  effect  of  three-dimensional  wave  fronts. 
The  determination  of  a  stiffness  profile  using  foe  SASW  method  consists  of 
three  phases:  data  measurement,  signal  analysis,  and  inversion,  each  of  which 
is  described  briefly  below. 
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Data  Measurement 


Signals  typically  are  recorded  from  two  sensors  as  shown  in  Figure  1  and 
stored  digit^ly.  Tlie  sensors  may  be  velocity  or  acceleration  transducers; 
acceleration  transducers  are  nom^ly  reserv^  for  cases  when  a  significant 
amount  of  high  frequency  energy  is  being  recorded.  If  the  source  is  impul¬ 
sive  (Figure  la),  a  time  window  with  the  full  length  of  the  transient  response 
is  collected  as  shown  in  Figure  2.  If  the  source  is  harmonic  (Figure  lb),  the 
cross-power  spectrum  phase  difference  and  coherence  between  the  source  and 
each  receiver  are  calculated  directly  using  a  dynamic  signal  analyzer,  averaged 
with  a  number  of  other  responses  using  the  same  configuration  and  source, 
and  stored  as  shown  in  Figure  3.  The  signal  for  the  source  is  taken  directly 
from  the  signal  generator  function  on  the  analyzer. 

The  most  common  procedure  for  measurement  is  to  keep  the  two  receivers 
separated  (in-line)  by  the  same  distance  that  the  first  receiver  is  from  the 
source  (common  midpoint  method)  as  depiaed  in  Figure  4.  Guidelines  for 
relationships  between  receiver  spacing  and  predominant  wavelengths  (function 
off)  have  been  proposed.  Heisey  (1982)  suggested  that  the  receiver  spacing, 
d  (=  dj  -  d,),  be  limited  to  the  range: 

^  ^  d  ^  2X  (2.2) 
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An  analytical  study  by  Sanchez-Salinero  (1987)  suggested  that  d2/d,  be  large 
(e.g.,  dj/d,  S  4).  However,  if  dj/d,  =  2,  then  the  wavelengths  considered 
should  be  limited  to: 

X  <  d,  (2.3) 


The  results  from  semi-analytical  studies  by  Sykora  (1993)  agree  with  the  work 
by  Sanchez-Salinero.  However,  from  a  practical  standpoint  for  surface  waves 
in  pavement  systems,  Rix  (1988)  suggests  that  the  range  of  acceptable  spac- 
ings  be  widened  to: 

^  ^  d  <  3X  (2.4) 


Signal  Analysis 

Signal  analysis  is  normally  performed  using  hardware  provided  with  the 
spectral  analyzer.  Data  for  harmonic  sources  are  immediately  converted  to 
cross-power  phase  and  coherence  and  averaged  as  described  in  the  previous 
section.  For  impulsive  sources,  the  time  records  of  displacement  are  first 
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Rgure  2.  Impulsive  signature  produced  by  hammer  striking  level  ground  at  receiver 
spacing  of  1 6  ft 
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Rgure  3.  Phase  of  cross*power  spectrum  and  magnitude  of  coherence  produced  by 
vibrator  on  level  ground  at  receiver  spacing  of  8  ft 
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Rgure  4. 


Arrangement  of  source  and  receiver  iitustrating  the  common  receiver  midpoint 
geometry  (Rix  1 988) 


averaged,  tli«i  the  data  analyzed  to  determine  the  variation  of  signal  phase  as 
a  function  of  source  frequency  and  develop  a  dispersion  curve  (variation  of 
phase  velocity  versus  wavelength). 

The  phase  velocity  and  wavelength  are  calculated  directly  from  the  real  and 
imaginary  parts  of  the  dynamic  displacements  (Nazarian  and  Stokoe  1985a). 

At  a  given  frequency,  /,  the  travel  time,  t,  between  two  receivers  is: 

t  =  (2-5) 

360  / 


where  and  ^  are  the  phases  at  points  1  and  2  (the  location  of  the  receivers) 
in  degrees  as  shown  in  Figure  5.  Then,  the  phase  velocity,  Q,  is: 

C  =  ^**^~*^j  =  360  f!^L^  (2.6) 

♦  t  (^2-^.) 

v^ere  the  wavelength,  X,  is: 
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X 


C* 

T 


(2.7) 


An  example  of  die  calculation  procedure  is  also  presented  in  Figure  S. 


Inversion 


Inversion  presently  is  an  iterative  process  whidi  is  indqiendent  of  wh^er 
die  data  were  collect^  using  an  inqiulsive  or  a  harmonic  source.  This  stqi 
involves  selecting  an  initial  profile  of  layer  thicknesses  and  moduli,  calculating 
the  dispersion  curve  using  matrix  solutions  (Haskell  1953  and  Thompson 
1950),  comparing  the  measured  and  matrix-calculated  dispersion  curves, 
modifying  Ae  profile  if  necessary,  then  repeating  the  process  until  a  close 
match  exists.  A  flowchart  of  this  process  is  shown  in  Figure  6.  The  inverse 
problem  does  not  necessarily  produce  a  unique  solution.  Automatic  matrix 
inversion  using  least  squares  is  used  at  UT  and  elsewhere  (Nazarian  and  Desai 
1993)  and  neural  networks  (Rix  and  Leipski  1991)  have  been  proposed  as  an 
alternative  to  the  original  matrix  scheme. 
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Figure  5.  Example  of  calculating  phase  velocity  from  unwrapped  phase 


Assume  a  Profile 


iKilliJ:- 


•  Shw  wave  Velocity 

•  Mass  Density 

•  Poisscm's  Ratio 


Final  Profile  is 
Assumed  to  Represmt 
the  Actual  Held  Conditions 
Accurately 


Figure  6.  Flow  chart  summarizing  the  inversion  procedure  (Rix  1 988) 
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3  Numerical  Approximation 
Method 


The  computer  code  vib3  is  a  research  tool  devdoped  at  WES  to  numeri¬ 
cally  approximate  surface  wave  propagation  in  complex  geosystems,  including 
sloping  ground.  This  code  was  aeded  using  principles  of  the  Finite  Elonent 
Mediod  (FEM)  and  written  using  the  Fortran  77  language  running  on 
U.S.  Army  and  Dqiartment  of  Defense  CRAY  supercomputers  at  WES 
(Sykora  and  Roesset  1992;  Sykora  1993).  It  has  been  validated  with  closed- 
form  solutions  for  axi-symmetric  problems.  Parametric  analysis  have  also 
been  conducted  to  define  ^propriate  ranges  for  system  and  problem  variables. 
Anodier  type  of  problem  examined  and  reported  dsewhere  (Sykora  1993)  is 
buried  layers,  bodi  continuous  and  discontinuous.  Salient  aspects  of  the  com¬ 
puter  code  are  presented  below. 


Formulation  and  Implementation 

The  primary  assumptions  used  to  develop  vib3  are  that  the  geometry  and 
boundary  conditions  of  the  system  and  the  distribution  of  material  properties 
are  planar  (2-D)  but  the  loads  can  be  non^lanar  (3-D).  This  set  of  conditions 
has  a  broader  range  of  applications  than  diat  for  plane  strain  (conunonly  used 
in  engineering  analysis)  while  circumventing  expensive  3-D  solution  methods. 
The  other  assumptions  are: 

a.  All  media  are  isotropic, 

b.  The  hysteretic  bdiavior  is  rq}resented  by  complex  moduli, 

c.  The  source  produces  vertical,  steady-state  excitation  on  the  ground 
surface  at  one  discrete  frequency,  and 

d.  The  distribution  of  loads  is  symmetric  in  the  out-of-plane  direction 
about  the  analysis  plane. 

Only  fixed  and  free  bpundary  conditions  were  used  for  this  study.  Fixed 
conditions  were  used  for  the  base  and  free  conditions  were  used  elsewhere. 
Boundaries  were  placed  sufficiently  far  from  the  source  and  points  of  interest 
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to  avoid  the  effects  of  any  reactions.  The  distances  to  boundaries  were 
sdected  based  on  parametric  analysis  of  boundary  locations. 

The  formulation  involves  two  primary  con^mnoits:  foe  condensation  of 
3-D  dynamic  stiffeess  matrices  to  equivalent  2-D  nuttrices  and  foe  rq)resenta' 
tkm  of  foe  distribution  of  loads  in  foe  out-of-plane  direction  using  a  Fourier 
expansion.  The  2-D  system  of  equations  are  first  solved  in  foe  frequency  and 
wave-number  domain;  inverse  Fourier  transforms  are  fora  performed  to  obtain 
foe  solution  as  a  function  of  out-of-plane  distance  and  time,  if  so  desired. 

The  FEM  is  a  numerical  analysis  technique  used  to  approximate  foe 
refuse  of  a  continuous  body  by  dividing  foe  domain  of  mterest  into  a  dis¬ 
crete  number  of  subdomains.  Boundary  conditions  and  extranal  forces  are 
imposed  at  discrete  nodes  where  foe  displacements  are  calculated.  Results  can 
be  interpolated  at  any  point  in  foe  body  through  foe  use  of  interpolation  func¬ 
tions.  In  general,  as  foe  subdomains  become  smaller,  foe  solution  converges 
to  that  of  foe  continuum. 

The  basic  equations  describing  a  steady-state  system  using  a  displacement 
formulation  are: 

,(K-«*M)U=P  (3.1) 


where 

K  s  "static"  stiffness  matrix 
(■)  s  circular  frequency  (rad/sec) 

M  =  lumped  mass  matrix 
y  =  vector  of  displacements 
P  =  vector  of  forces 

Calling  S  foe  dynamic  stiffeess  matrix  defined  by: 

S  =  K  -  M  (3*2) 


fora 


S  U  =  P 


(3.3) 


Modification  to  constitutive  model 

Soil  is  an  indastic  material>-raergy  dissipates  from  friction  braveen  parti¬ 
cles  as  waves  travel  through  U.  This  loss  of  energy  is  rq)resented  through  foe 
use  of  material  damping.  Difrrarat  mathematical  modds  are  used  to  ^proxi¬ 
mate  damping  in  foe  governing  equations.  One  form,  called  hysteretic  damp¬ 
ing,  is  frequency  indqpendent.  This  form  of  damping  can  be  introduced  into 
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the  fonnuladon  for  frequency-domain  analyses  tfarou^  die  Correspondence 
Principle  (Wolf  1985).  This  principle  states  that  the  elastic  stifim^  (in  this 
case  shear  modulus,  G)  is  rqilaced  by  a  conqilex  stiffoess  to  otoun  the 
damped  solution.  The  following  relationship  is  commonly  used  to  model 
linear-hysteretic  behavior  for  small  shear  strains  (and  small  values  of 
damping): 

G  •  *  G  (1  ♦  2i(S) 


where 


G*  =  complmc  shear  modulus 
fi  =  damping  ratio 

The  results  of  this  study  are  expected  to  be  applied  at  distances  greater  than 
one  wavelength  from  the  source  where  shear  strains  are  small  for  synthetic 
sources.  The  magnitude  of  damping  is  considered  to  be  independent  of  strain 
(Hardin  and  Dmevich  1972;  Jolmston,  Toksoz,  and  Timur  1979;  and  Toksoz, 
Johnston,  Timur  1979)  for  the  levels  of  shear  strains  expected. 


Fourier  superposition 

Fourier  superposition  is  a  three-stqp  solution  process  for  linear  systems  that 
involves  a  forward  transformation  into  a  wavenumber  domain  (either  spatial  or 
temporal),  the  calculation  of  a  solution  vector  of  displacements  (Eq.  3.3)  at  a 
number  of  incrments,  and  the  determination  of  the  total  solution  through  an 
inverse  transformation  of  all  incremratal  solutions.  Fourier  superposition 
sq)plied  in  both  the  time  and  y-spatial  domains  leads  to  the  transform  pair: 

p(in,tf)  “1*1*  dt  dy  (3-5) 

iKy..)  -  ±  I 


I  **  p(m,(i))  dco  dm 


(3.6) 


where 

p  -  distribution  of  forces  with  respect  to  y-direction  and  time 
p  =  distribution  of  forces  in  wavenumber  domain 
m  =  wavenumber  (spatial  circular  frequency)  in  y-direction 

Making  the  load  vector  specific  to  steady-state  vibrations  with  constant  ampli¬ 
tude,  dien: 
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(3.7) 


P(W|).  -  J/J  P(y)  e*^  dy 


^1 

“  ir  1--  e-*^  dm 

for  a  specific  u.  The  corresponding  equations  for  displacements  are: 

u(m)^  *  |.r  “(y) 


u(y.t) 


^  L.  e-J-  dm 


(3.8) 


(3.9) 

(3.10) 


where 

u  =  distribution  of  displacements  with  respect  to  y-direction 
u  =  distribution  of  displacements  in  wavenumber  domain 


Element  condensation 

The  process  of  element  condensation  is  the  key  aspect  of  the  reduction  of 
computational  effort.  Element  condensation  refers  to  the  process  of  reducing 
the  number  of  degrees  of  freedom  by  relating  points  adjacent  in  the 
y-direction  using  the  functional  relationship  of  the  Fourier  expansion.  The 
dq)endent  degrees  of  freedom  are  then  eliminated  by  expressing  them  in  terms 
of  the  in-plane  degrees  of  freedom. 

Finite  element  condensation  has  been  used  successfully  by  a  number  of 
previous  investigators.  This  strategy  was  explicitly  proposed  for  axi- 
symmetric  problems  by  Winnicki  and  Zienkiewicz  (1979)  and  Lai  and  Booker 
(1991)  and  for  3-D  formulations  by  Runesson  and  Booker  (1982,  1983)  and 
Lin  and  Tassoulas  (1987).  Element  condensation  was  reported  by  Kang 
(1990)  and  Hanazato  et  ai.  (1991)  for  wave  propagation  studies  in  horizontally 
layered  pavement  systems.  For  this  study,  the  degrees  of  freedom  corre¬ 
sponding  to  the  nodes  outside  of  the  x-z  plane  are  eliminated.  Each  node  in 
the  two-dimensional  mesh  maintains  three  degrees  of  freedom. 

Consider  an  arbitrary  discretized  model  of  a  slope  system  that  meets  the 
requirement  of  uniform  geometry  and  material  properties  in  one  direction  such 
as  that  shown  in  Figure  7.  The  coordinate  system  is  chosen  to  have  the 
z-direaion  positive  down  and  the  other  in-plane  direction  to  be  x.  (Consider 
diree  verti(^  planes  sq>arated  by  a  distance  of  Ay  at  some  arbitrary  location 
along  the  geosystem.  The  dynamic  stiffness  matrix  for  a  single  element 
betwera  two  slices  can  be  partitioned  as: 
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Figure  7.  Example  of  a  physical  system  for  purposes  of  extracting  a  slice  of  finite 
elements 
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(3.11) 


where  the  subscripts  *  1"  and  "2"  refer  to  the  dq^rees  of  freedom  on  the  posi¬ 
tive  and  n^ative  face  in  die  y-direction,  reflectively.  The  assemblage  of  the 
dynamic  equations  for  any  two  finite  elements  adjacent  in  the  y-direction,  as 
shown  in  Figure  8b,  can  be  reduced  by  canceling  the  time-de(Mendent  exponen¬ 
tial  term  on  each  side  to: 

Sn  &i2 
$21  Sa+S,, 

0  821 

where 

y  =  2-D  vector  of  displacements  on  designated  face 
P  =  2-D  vector  of  forces  acting  on  designated  face 
"  +  "  denotes  element  in  -t-  y-direction  (from  Ay  to  0) 
denotes  element  in  -  y-direction  (from  0  to  -Ay) 

”a"  denotes  degrees  of  freedom  on  face  a  (i.e.,  at  y  =  0) 

"b*  denotes  degrees  of  freedom  on  face  b  (i.e.,  at  y  =  +Ay) 

"c"  denotes  degrees  of  freedom  on  face  c  (i.e.,  at  y  =  -Ay) 

Using  the  Fourier  expansion  described  earlier,  forces  and  displacements  are 
expressed  as: 


Pirn)  * 

1  ^  P(y)  e*^  dy 

(3.13) 

C  (w)  » 

1  *  U(y)  e*^  dy 

(3.14) 

where  P  and  C)  are  used  to  represent  vectors  of  nodal  forces  and  displace¬ 
ments,  respectively,  in  m  space.  Rewriting  Equation  3.12  to  incorporate  the 
Fourier  expansion  of  loads: 


fr^(m) 

^,(m) 

^.(m) 


(3.15) 
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In  die  transform  (m)  space,  die  displacmnents  on  the  *b"  and  *c*  faces  are 
related  to  die  di^lacements  on  die  *a”  dice  at  any  instant  in  time  by  the  sim¬ 
ple  relationships: 

tfc  (m)  -  U,(m)  e-^  (3.16) 

Cr,  (m)  =  U,(m)  e^y  (3- 17) 

Defining: 

8(m)  »  s^,  e‘*^>'  ■*■  (  *11  +  S22  )  ♦  i,2  (3-18) 


which  is  not  symmetric— the  stiffness  terms  corresponding  to  displacements  in 
the  y-dire(^ion  are  out  of  phase.  The  stiffness  matrix  is  made  symmetric  by 
multiplying  the  rows  and  dividing  the  columns  corresponding  to 
y-displacements  by  i  and  solving  in  terms  of  iv  where  v  is  the  y-displacement. 
Equations  3.16  and  3.17  can  be  substituted  into  Equation  3.15  to  get  the 
system  of  equations  for  the  equivalent  two-dimensional  system  shown  in 
Figure  8c: 

S(m)  0,(m)  =  ^,(m)  (3.19) 


This  formulation,  then,  allows  the  three-dimensional  element  with  a  two- 
dimensional  geometry  to  be  represented  with  an  equivalent  2-D  element. 


Validation  and  Parametric  Analysis 

Validation  studies  were  conducted  to  prove  that  the  formulation  and  com¬ 
putet  implementation  are  sound  for  the  limited  class  of  problems  considered. 
The  general  accuracy  of  calculated  displacements  in  terms  of  amplitude  and 
phase  was  evaluated.  A  brief  presentation  of  the  results  of  the  validation  and 
parametric  analysis  are  described  below. 

The  best  form  of  validation  consists  of  comparisons  with  exact  mathemati¬ 
cal  relationships  for  a  suite  of  problems  representing  a  wide  spectrum  of 
possible  in-situ  conditions.  Comparisons  with  measured  data  from  full-scale 
or  prototype  testing  provide  a  constructive  means  to  confirm  findings  when 
co^ucted  undo*  certain  controlled  conditions.  Comparisons  with  measured 
data  are  not  iqipropriate  as  the  primary  and  only  means  of  validation, 
howevCT.  Comparisons  with  ofoer  numerical  approximations  are  even  less 
appropriate  for  validation.  Validation  of  wb3  through  comparisons  with  ana¬ 
lytical  results  is  possible  only  for  the  simplest  class  of  planar  geometry— a 
horizontally  layered  system  extending  to  infinity.  Green’s  function  solutions 
formulated  for  axi-symmetric  problems  (with  disk  loads)  by  Kausel  (1981) 
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were  used  exclusively  with  point  and  square  loads  (in  plan).  The  same  total 
area  and  a  total  load  of  unity  were  used  to  adequately  equate  problems  with 
different  load  configurations. 

The  dynamic  vertical  displacements  along  the  surface  and  sufficiently  away 
from  the  load  are  of  primary  interest.  Vertical  vibrations  are  normally  mea¬ 
sured  in  non-destructive  testing  techniques  such  as  the  SASW  method.  Sen¬ 
sors  are  easy  to  align  in  the  vertical  position  regardless  of  the  surface  slope. 
Measurements  are  likely  to  be  made  both  perpendicular  and  parallel  to  the 
structure  of  the  system  (x-  and  y-directions,  respectively).  Ifierefore,  the 
results  are  presented  in  terms  of  the  variations  of  real  and  imaginary  compo¬ 
nents  of  dynamic  displacement  in  both  horizontal  directions— the  x-direction 
(calculated  at  nodes  on  the  y  -  0  plane)  and  the  y-direction  (calculated  at  the 
node  nearest  the  centroid  of  the  load  and  expanded  out  in  the  y-direaion). 

The  variation  of  phase  with  distance  is  added  to  provide  a  means  for  closer 
examination. 

Comparisons  between  the  numerical  approximations  and  closed-form  solu¬ 
tions  at  distances  from  the  source  between  O.S  X  and  S.O  X  are  excellent  for 
nearly  all  of  the  cases  considered.  Examples  of  the  close  comparison  are 
shown  in  Figures  9  through  11.  The  vertical  components  of  dynamic  dis¬ 
placement  at  the  surface  and  at  a  depth  of  0.8  X  are  shown  in  Figures  9  and 
10;  respectively.  Although  not  considered  further  in  this  report,  the  horizon¬ 
tal  components  of  displacement  at  the  surface  are  shown  in  Figure  11.  The 
data  in  these  figures  as  well  as  many  other  comparisons  not  shown  confirm 
that  the  formulation  and  computer  implementation  are  accurate. 

Parametric  analyses  were  conducted  to  assess  the  sensitivity  of  the  formula¬ 
tion  and  numerical  solution  to  system  variables  and  to  recommend  thresholds 
for  these  variables.  The  formulation  was  found  to  be  sensitive  to  the  follow¬ 
ing  system  parameters:  the  spatial  increments  of  discretization  in  all  three 
directions  and  the  total  length  of  discretized  space  in  the  y-direction.  The 
results  are  also  moderately  sensitive  to  reasonable  ranges  of  Poisson’s  ntio 
and  damping  for  soils.  The  solution  accuracies  are  not  sensitive  to  the  width 
of  the  load  using  reasonable  bounds  of  load  width,  soil  velocities,  and  operat¬ 
ing  frequencies.  The  calculated  displacements  were  shown  to  be  slightly  more 
accurate  as  damping  ratio  was  increased  from  2  to  S  percent. 

Comparisons  were  also  made  between  the  time  and  memory  requirements 
to  execute  vib3  and  an  equivalent  solution  model  using  a  commercial  3-D 
FEM  software  package  {ABAQUS).  The  semi-analytical  method  was  found  to 
generally  be  25  times  faster  and  use  only  one-fourth  the  CPU  memory  space. 
Therefore,  the  new  code  is  an  accurate  and  economical  alternative  to  expen¬ 
sive  3-D  modeling. 
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Model  1  Ax/X  =  0.10  Ay/A  =  0.05 

^  =  0.05  /  =  3.0  Hz  YTOT/X  =  ±  6.4 


Figure  9.  Comparison  between  calculated  vertical  displacements  and  Green's  function 
solution  along  surface 
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VERTICAL  DYNAMIC  DISPLACEMENTS 


Model  1  Ax/A  =  0.20  Ay/A  =  0.05 

P  =  0.05  /  =  3.0  Hz  YTOT/A  =  ±  6.4 

Depth/A  =  0.8  NM  =  256  XLOAD/A  =  ±  0.01 


o 


HORIZONTAL  DISTANCE  (in  X  ) 


Figure  10.  Comparison  between  calculated  vertical  displacements  and  Green's  function 
solution  at  depth  of  0.8  >1 
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Figure  1 1 .  Comparison  between  calculated  horizontal  displacements  and  Green's  function 
solution  along  surface 
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4  Results  of  Computations 


This  chapter  is  used  to  describe  the  results  of  analyses  compiled  with  the 
semi-analytical  method.  The  variations  of  dynamic  vertical  displacements  in 
the  x-direction  (direction  of  maximum  topographic  gradient)  ai^  the 
y-direction  (direction  of  zero  topographic  gradient)  e  calculated  and  ana¬ 
lyzed  to  define  trends  that  can  ^  used  with  inversion  schemes  to  derive  mate¬ 
rial  stiffness.  The  slope  was  assumed  to  extend  to  infinity. 

Several  researchers  have  examined  die  effect  of  surface  discontinuities  on 
die  amounts  of  surface  wave  energy  transmitted  and  refleaed  using  laboratory 
measurements.  The  collection  of  Aeir  results  is  somewhat  inconsistent  but 
dott  show  a  strong  dependence  on  the  slope  angle.  A  few  have  also  exam¬ 
ined  the  effect  of  surface  discontinuities  on  phase.  These  results  indicate  that 
die  effect  of  slope  angle  can  be  easily  avoided  or  taken  into  account. 

The  variables  defining  the  problem  are  shown  in  Figure  12  with  the  slope 
angle,  a,  and  the  load  s^-back  distance,  x^;,  being  the  primary  variables; 
infinite  slope  conditions  were  modeled.  Square  loads  were  us^  throughout 
and  can  be  rdated  to  circular  load  areas  by  use  of  equivalent  total  loads. 
Calculations  were  made  using  six  common  engineering  slope  ratios  (horizon¬ 
tal :  vertical)— 1:1,  1.5:1,  2:1,  2.5:1,  3:1,  and  6:1.  The  load  was  placed  at  the 
top  and  base  of  the  slope  widi  these  ratios  to  produce  a  =  135,  146,  153, 

158,  162,  171,  189,  198,  202,  207,  214,  and  225  degrees.  (An  a  of  180 
degrees  corresponds  to  the  level  ground  conditions  considered  in  previous 
chsq>ters.)  The  center  of  the  load  was  used  as  the  reference  point  for  distance 
measuremmts.  Unit-less  dimensions  are  used  for  convenience.  The  shear 
wave  velocity  was  fixed  at  1000  with  5  percent  danq)ing.  The  mass  density 
and  Poisson’s  ratio  were  4.00  and  0.40,  respectively,  producing  a  equal 
to  940. 

The  results  of  computations  were  initially  evaluated  to  determine  the  effect 
of  system  parameters.  Thresholds  defined  by  Sykora  (1993)  for  level  ground 
conditions  were  generally  found  to  be  ^propriate  for  sloping  ground  condi¬ 
tions.  Free  and  fixed  (rigid)  boundaries  of  foe  discretized  FEM  models  were 
placed  at  sufficient  distances  to  avoid  any  effects  in  foe  regions  of  comparison 


One  initial  finding  was  that  foe  phase  velocities  calculated  from  measure- 
mrats  along  an  infinite,  homogeneous  slope  are  independent  of  foe  harmonic 
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XLOAD  e 


Figure  12.  Schematic  of  problem  variabies  used  to  analyze  surface  wave  propagation 
sloping  ground 


frequency  of  excitation  (i.e.,  non-dispersive).  This  finding  is  consistent  with 
levd  ground  conditions  and  greatly  simplifies  the  analysis  of  data  for  sloping 
ground  conditions. 

The  variations  of  vertical  displacement  and  phase  at  a  slope,  for  a  slope 
ratio  of  1.5:1  (a  =  146  d^rees  the  load  is  at  the  crest  of  the  slope  aid 
a  =  214  d^ees  ^en  the  load  is  at  the  base  of  the  slope)  are  shown  in 
Figures  13  and  14,  respectively.  The  relative  location  of  the  discontinuity 
along  die  abscissa  is  shown  as  a  vertical  dashed  line.  The  level-ground  solu¬ 
tion  (Kausd  1989)  is  dso  shown  for  comparison. 

The  primary  effect  of  die  slope  discontinuity  is  a  phase  shift  for  waves 
travding  in  die  direction  of  maximum  topographic  gradient  that  occurs 
between  die  discontinuity  and  a  distance  of  3X/4  beyond  (when  x^:  =  1.0  X). 
The  phase  diffn’ence  becomes  essentially  constant  at  greater  distances.  The 
phase  has  a  marked  lag  from  the  levd  ground  solution  at  distances  beyond  the 
discontinuity  idien  a  =  146  degrees  and  a  marked  phase  lead  when 
a  =  214  d^ees.  The  calculated  solutions  in  die  y-direction  along  die  crest 
and  base  correspond  wdl  to  die  levd-ground  solution,  regardless  of  a. 

The  amount  of  phase  shift  at  these  two  slope  angles  was  examined  by 
comparing  die  total  (unwrapped)  phase  widi  the  total  phase  from  levd  ground 
(Green’s  function)  solutions,  and  the  total  idiase  produced  by  a  uniform 
increase  of  360  d^ees  per  wavdength  and  passing  forougb  the  origin, 

These  differences  are  shown  in  Figures  IS  and  16,  respectivdy.  The  is 
used  to  calculate  di^iersion  curves  for  the  SASW  method  and  a  comparison 
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SURFACE  DISTANCE  (in  X  ) 


Rgure  14.  Variation  of  displacements  and  phase  for  o  ■  214  degrees 
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a  =  146*  X,,  /X  =  1.00  XLOAD/X  =  ±  0.01 

^  =  0.05  /  =  3.0  H2 


SURFACE  DISTANCE  (in  A  ) 

a.  a  «  146  degrees 


a  =  214*  /X  =  1.00  XLOAD/X  =  ±  0.01 

^  =  0.05  /  =  3.0  H2 


b.  a  214  degrees 


Figure  1 5.  Difference  in  phase  between  sloping  ground  and  level  ground  conditions 
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DIFFERENCE  IN  PHASE,  n  DIFFERENCE  IN  PHASE, 


a  =  146  /X  =  1.00  XLOAD/X  =  ±  0.01 

0  =  0.05  /  =  3.0  Hz 


SURFACE  DISTANCE  (in  X  ) 


1 46  degrees 


a  =  214*  x,^/X  =  1.00  XLOAD/X  =  ±  0.01 

^  =  0.05  /  =  3.0  Hz 


SURFACE  DISTANCE  (in  X  ) 


|b.  g  =  214  degrees _ 

Rgure  1 6.  Difference  in  phase  between  sloping  ground  and  a  uniform  increase  of 
360  degrees  per  wavelength 
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between  ^  and  ^  is  shown  in  Figure  17.  The  phase  diffidence  between  the 
two  vertical  values  s^proaches  60  d^ees  at  distances  less  than  SX  when 
»  =  0.40. 

The  difference  in  phase  between  sloping  and  level  ground  conditions  in 
Figure  15  shows  the  relatively-small  nuignitude  of  phase  shift.  The  variation 
of  phase  difference  after  the  ^ift,  beyond  distances  of  about  3X/4  from  the 
discontinuity,  is  relatively  constant  when  a  =  146  degrees.  The  small  fluctu¬ 
ations  that  exist  are  similar  to  those  for  the  solution  along  the  crest.  The 
overall  variation  of  phase  in  die  x-direction  when  a  =  214  degrees  is  much 
larger  and  less  consistent.  Along  the  base  and  crest  there  is  no  noticeable 
phase  difference  and  the  magnitude  of  fluctuations  is  small. 

The  difference  in  phase  between  sloping  conditions  and  ^  is  shown  in 
Figure  16.  Both  sets  of  data  have  shifted  down  relative  to  the  data  shown  in 
Figure  IS  because  of  the  difference  bttween  and  4m  noted  previously. 

The  differences  along  the  crest  and  down  the  slope  for  or  =  146  degrees  again 
are  both  relatively  constant.  These  data  strongly  suggest  that  measurements 
can  be  reliably  used  in  both  the  x-direction  (down  the  slope)  and  the 
y-direction  to  determine  soil  stiffness  using  the  SASW  method,  especially  if 
the  receivers  are  placed  at  a  distance  greater  than  3X/4  beyond  the  discontinu¬ 
ity  to  avoid  the  phase  shift. 

The  results  when  a  =  214  degrees  are  not  as  consistent.  The  phase  differ¬ 
ence  i^pears  to  vary  considerably  at  surface  distances  beyond  the  initial  phase 
shift.  The  general  trend  at  distances  greater  than  2.S  X  is  for  the  phase  to 
decrease  uniformly  (lag)  as  a  ratio  of  distance  until  a  sudden  shift  occurs 
(e.g.,  at  3.4  X).  This  type  of  response  is  not  desirable  for  methods  such  as 
the  SASW. 

The  variation  of  phase  from  a  larger  perspective  indicates  that  the  differ¬ 
ences  shown  in  Figures  IS  and  16  may  not  be  of  importance  from  a  practical 
standpoint.  The  variation  of  unwrapped  phase  with  distance  for  most  slope 
angles  (minor  discontinuities  not  included)  are  shown  in  Figures  18  and  19. 
Also  shown  in  Figure  18  are  lines  corresponding  to  a  least  squares  linear  fit  of 
data  between  surface  distances  of  1.0  and  5.0  X  for  each  slope  angle. 

The  data  in  Figure  18  shows  that  the  distribution  of  unwrtqiped  phase  is 
fairly  consistent  and  unique  when  a  <  180  degrees.  The  best-fit  segments  at 
these  angles  effectively  represent  a  variation  of  phase  between  360  and 
350  degrees  per  wavelength.  The  increased  variation  in  phase  when 
a  >  180  degrees  is  also  seen  in  Figure  18  although  the  best-fit  linear  seg¬ 
ments  are  very  similar  and  are  close  to  4^.  The  error  associated  with  this 
variation  is  small  and  can  be  minimized  by  using  large  spacings.  The  data  in 
Figure  19  clearly  shows  the  independence  of  phase  change  along  the  crest  and 
base  to  a  and  how  close  the  phase  is  to  43(0  (shown  as  a  line  segment  in 
Figure  19). 

The  phase  shifts  for  all  a  considered  using  the  differences  shown  in 
Figure  16  were  plotted  in  Figure  20.  The  difference  in  phase  is  not  zero  at 
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PHASE,  ^  (degrees) 


Rgure  17.  Difference  in  Green's  function  solutions  for  horizontal  and  vertical  components 
(Kausel  1989) 
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PHASE,  ^  (degrees) 


Figure  1 8.  Variation  of  phase  in  direction  of  greatest  topographic  gradient 
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PHASl,  0  (degrees) 


Rgure  19.  Variation  in  phase  in  direction  of  zero  topographic  gradient 
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Figure  20.  Variation  of  phase  shift  with  slope  discontinuity  angle 


a  =  180  degrees  which  reflects  the  differrace  between  0^  and  These 
results  are  smaller  than  "first-order  theory”  (no  surface  wave  reflections 
included)  reported  by  Hudson  and  Knopoff  (l%4).  Although  it  appears  possi¬ 
ble  to  predict  the  phase  shift  caused  by  a  surface  discontinuity,  calculated 
phase  velocities  are  likely  to  be  more  accurate  if  this  phase  shift  is  avoided  in 
the  analysis.  The  difference  in  phase  as  a  function  of  a  was  found  to  be 
independent  of  x^;  as  long  as  Xi£  >  0.5  X. 

The  real  and  imaginary  parts  of  peak  displacements  at  distances  beyond  the 
discontinuity  amplify  as  a  decreases  from  180  degrees  and  deamplify  as  a 
increases  above  180  degrees.  The  ratio  of  peak  vertical  displacements  as  a 
function  of  a  to  peak  vertical  displacements  for  level  ground  are  shown  in 
Figure  21.  This  ratio  is  a  function  of  x^;  and  tends  to  be  a  maximum  when 
Xix;  is  about  1.0  X.  A  standard  set-back  distance  of  1.0  X  was  used  for  most 
analyses. 

The  combined  results  suggest  that  measurements  made  along  the  crest  or 
base  of  a  slope  can  safely  be  processed  using  existing  SASW  inversion 
schemes  to  obtain  stiffness  profiles.  l%ase  velocities  rq)resenting  measure¬ 
ments  down  a  slope  can  also  be  reliably  obtained  from  vertical  measurements. 
The  best  method  would  be  to  place  receivers  beyond  3X/4  thus  avoiding  any 
interpretation.  Much  more  variation  is  expeaed  for  measurements  made  up  a 
slope.  For  homogeneous  systems,  the  phase  velocity  will  closely  !q)proximate 
Vr. 


34 


Chapter  4  Reaults  of  Computations 


Figure  21 .  Variation  of  ratio  of  peak  vertical  displacements  with  slope  discontinuity 
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Field  Measurements 


Field  measurements  of  vertical  dynamic  displacements  on  level  ground  and 
sloping  ground  were  made  on  die  Reservation  in  Vicksburg,  Mississippi, 
during  the  month  of  August,  1993.  These  data  allow  comparisons  between 
the  results  of  numerical  analyses  and  physical  measurements,  particularly  since 
exact  solutions  do  not  exist  for  sloping  ground.  A  secondary  purpose  for 
conducting  the  field  program  was  to  waluate  the  effect  of  seismic  sources  and 
sensor  spacings  on  measured  displacements.  Discussions  about  the  response 
produced  by  different  combinations  of  these  two  factors  are  included  bdow. 

The  site  is  located  about  100  yards  northeast  of  the  Casagrande  Building 
and  adjacent  to  Porter’s  Chapel  Road.  It  consists  of  a  road  cut  with  a  maxi¬ 
mum  height  of  25  ft  and  slope  angles  (using  the  convention  presented  in  Fig¬ 
ure  12)  between  about  155  and  165  degrees  and  a  concrete  bridge  with  a  steel 
support  system.  A  site  map  is  provided  in  Figure  22  and  a  photograqph  of  the 
bridge  is  shown  in  Figure  23.  The  southern  portion  of  the  site  is  flat  and 
level  but  a  short  road  embankment  ^ich  serves  as  an  ^proach  for  the  bridge 
exists  on  the  eastern  margin.  The  crest  of  the  cut  slope  lies  in  an  ESE-WNW 
direction. 

Natural  soil  consists  of  Pleistocene-age  loess  indigenous  to  the  Vicksburg 
area.  This  loess  is  uniformly  graded  and  typically  has  a  low  unit  weight. 
Groundwater  does  not  exist  within  the  d^tbs  tested.  During  testing,  the 
ground  surface  was  very  resistant  to  penetration  which  is  attributed  to 
desiccation. 

Some  shallow  fill  may  also  exist  across  the  site,  the  consequence  of  site 
grading  during  construction.  An  embankment  fill  was  used  beneath  tiie 
county  road  as  an  ^proach  to  the  bridge. 


General  Procedure 

A  field  testing  program  was  designed  to  determine  the  stiffness  profile  on 
levd  ground  (free  fidd)  to  a  dq)th  of  about  25  ft  and  to  measure  dynamic 
vertical  di^lacements  along  die  crest  and  down  the  slope  at  two  locations. 
The  general  locations  of  the  arrays  are  shown  in  Figure  22.  Both  inq)ulsive 
and  harmonic  wave  sources  were  used.  Impulsive  sources  were  hammer 
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Figure  22.  Site  map  of  WES  test  site 
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Figure  23.  View  of  bridge  at  WES  test  site  looking  east  with  area  of  seismic  measurements  on  sloping  ground  in  foreground 


blows,  dropping  of  a  boulder,  and  activmioo  of  a  Dynatest  Falling  Weight 
Deflectometer  (FWD).  Greater  energy  is  typically  required  as  the  spacing 
increases  as  was  die  case  at  diis  site.  The  harmonic  excitation  was  produced 
using  a  SO-lb  Goodman  electro-magiietic  vibrator  with  a  platm  diameter  of 
6  in.  A  Hewlett  Packard  3S62A  dynamic  signal  analyzer  was  used  to  produce 
die  uqiut  signal  and  to  collect,  analyze,  and  store  sig^s.  An  overall  view  of 
the  equipment  used  for  nwasurements  with  harmonic  loads  is  shown  in 
Figure  24. 

Averaging  of  signals  was  used  for  collecting  data  from  both  source  types. 
For  impulsive  sources,  five  time-domain  signals  were  averaged  and  stored 
digitally  using  die  HP  spectrum  analyzer  (example  shown  previously  in 
Figure  2).  For  harmonic  source  measurements,  the  vibrator  was  programmed 
to  run  a  swept-sine  test  between  frequencies  of  25  and  325  Hz  at  increments 
of  0.375  Hz.  Five  cross-correlations  and  coherency  relationships  between  the 
source  and  each  receiver  were  averaged  at  eadi  frequency  and  stored  (example 
shown  previously  in  Figure  3). 

Dispersion  curves  were  developed  from  data  respective  of  each  location 
and  direction.  Inversion  techniques  devdoped  by  Nazarian  (1984)  and 
adapted  by  WES  (Alexander  19^)  w^e  employed  to  develop  dispersion 
curves  for  data  collected  with  impulsive  sources  (refn*  to  Chapter  2).  The 
variation  of  phase  and  coherence  calculated  from  a  cross  power  spectrum 
between  die  two  averaged  signals  with  frequency  are  the  primary  relationships 
used  for  this  process.  Dispersion  curves  for  data  collected  using  the  harmonic 
source  were  easily  derived  by  subtraaing  values  of  (unwrapped)  phase 
between  the  two  receivers  and  using  Equation  2.6. 


Level  Ground  Measurements 

The  testing  configuration  for  level  ground  measurements  is  shown  in  Fig¬ 
ure  25  and  follows  the  common-midpoint  procedure  with  dj/d,  =  2  (except 
for  Test  FIO  where  dj/d,  =  1.7).  Velocity  transducers  (Mark  Products  L-22) 
with  a  resonant  frequency  of  2  Hz  were  positioned  in  an  attempt  to  avoid  the 
effects  of  the  nearby  slope,  pavements,  structures.  Both  impulsive  and 
harmonic  sources  were  used  to  measure  surface  wave  displacements  in  level 
ground.  However,  the  dispersion  curve  and  stiffness  profile  were  derived 
only  from  measurements  using  the  impulsive  source. 

The  presence  of  an  embankment,  even  a  short  embankment,  was  found  to 
introduce  some  scattering  of  the  signal.  Consequently,  waves  traveling  east 
(Tests  FI,  F3,  and  F6)  were  not  used  to  evaluate  stiffness.  The  variation  of 
phase  and  cohwence  with  frequency  calculated  from  hanuner  blows  at  a  spac¬ 
ing  of  8  ft  in  opposite  directions  (Jests  F3  and  F4)  are  shown  in  Figure  26. 
The  data  from  Test  F3  correspond  to  sensors  being  located  between  the  source 
and  the  county  road  embankment;  the  data  from  Test  F4  correspond  to  sensors 
being  located  in  the  direction  away  from  the  road  and  source.  The  sources 
were  60  ft  and  36  ft  from  the  base  of  the  embankment  for  tests  F3  and  F4, 
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Rgure  25.  Testing  configuration  for  level  ground  measurements 


Cross  Power  Spectrum  Cross  Power  Spectrum 
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e  26.  Phase  difference  and  coherence  between  receivers  at  8-ft  spacing  produced  by  averaged  hammer  blows  in  level  ground 


rcqMCtivdy.  The  |riiase  rdMiowhy  for  wavm  triveiim  away  ftom  the 
ftmhmnicmfmt  is  »«■«*  finopAer  and  die  oohercoce  is  generally  hi^er  for 
frequencies  between  25  and  12S  Hz.  SimUar  findings  were  made  at  otfaff 
HMcings. 

A  significant  difference  was  found  to  exist  in  the  energy  imparted  by  the 
boulder  and  the  FWD  for  die  level  ground  measurements.  Corresponding 
signals  for  Tests  FIO  and  FI  1  are  shown  in  Figure  27.  The  coherence  pro¬ 
duced  using  the  FWD  at  frequencies  greater  than  IS  Hz  are  quite  large,  much 
bett«  than  the  coherence  produced  using  the  boulder.  A  conqiarison  of  phase 
relationships  is  less  conclusive,  however.  Based  on  this  comparison,  it  would 
aiqiear  diat  the  FWD  is  a  good  source  for  large  spacings  on  level  ground, 
mudi  better  than  a  boulder. 

The  level  ground  site  was  found  to  have  a  non-uniform  profile  of  stifhiess 
(which  was  the  assumption  for  numerical  calculations).  The  dispersion  curve 
for  level  ground  measurements  using  all  three  impulsive  sources  is  shown  in 
Figure  28  widi  the  individual  test  data  indicated.  A  best-fit  relationship  used 
for  inversion  is  also  shown.  The  phase  velocity  generally  increases  as  a  func¬ 
tion  of  dqith  to  the  power  of  0.27.  A  comparison  between  final  iterated 
points  along  the  dispersion  curve  and  measured  data  is  made  in  Figure  29. 

The  stifih^  profiles  produced  using  the  best-fit  relationship  of  the  dispersion 
curve  assummg  a  Poisson’s  ratio  of  0.40  and  a  density  of  100  pcf  is  shown  in 
Figure  30.  The  increments  in  modulus  with  dq>th  produce  a  relatively  uni¬ 
form  increase  with  dqtth  (2,200  psf  per  ft  dq)th)  excq)t  between  depths  of  6 
to  1 1  ft  where  the  increase  is  sli^tly  less. 


Sloping  Ground  Measurements 

The  locations  of  the  slope  arrays  were  shown  previously  in  Figure  22.  The 
testing  configuration  for  slope  measurements  and  additiomd  photographs  of 
diese  areas  are  shown  in  Figures  31  dirough  33.  The  slope  angles  measured 
using  a  level  survey  at  Locations  A  and  B  are  160.5  and  163  degrees,  respec¬ 
tively.  The  testing  configuration  was  selected  to  avoid  measuring  phase  shifts 
near  the  discontinuity  and  allow  the  use  of  a  wide  range  of  frequencies. 

Twdve  velocity  tran^ucers  (Mark  Products  L20  DX’s)  with  a  resonant  fre- 
qumicy  of  14  Hz  were  positioned  vertically  at  locations  down  the  slope  and 
^ong  die  crest  (toward  the  other  array  location)  at  1-ft  spacings  between 
distances  of  10  and  21  ft  from  the  center  of  the  source  platen.  (At  site  A, 
only  10  geophones  were  used  down  the  slope.  In  addition,  the  data  at  a  sur¬ 
face  distance  of  16  ft  were  inadvotendy  not  stored.) 

In  most  instances,  the  phase  and  coherence  for  measurements  using  a  har¬ 
monic  source  for  sloping  ground  conditions  are  nearly  ideal  for  processing. 
One  such  set  of  phase  and  coherence  data  is  shown  in  Figure  34.  The  phase 
is  nearly  linear  and  the  coherence  is  usually  greater  than  0.99  excqit  at  fre¬ 
quencies  less  than  50  Hz  and  greata*  than  250  Hz.  Aliasing  of  phase  relation- 
^ips  typically  began  to  occur  at  wavelengdis  less  than  3  ft  (frequencies 
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VELOCITY,  C  (fps) 
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Figure  28.  Dispersion  curve  derived  from  level  ground  measurements 


e  29.  Final  iterated  dispersion  curve  in  comparison  to  best  fit  of  level  ground  measurements 


PHASE  VELOCITY,  C  *  (fps) 
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Chapter  5  Raid  Maaeuremanta 


49 


Figure  32.  View  of  test  configuration  for  measurements  made  downslope  looking  southeast 


Figure  33.  View  of  test  configuration  for  measurements  made  along  the  crest  looking  northwest 


difference  and  coherence  between  source  and  receiver  for  slope  measurements  and  harmonic  source 


greater  then  about  125  Hz).  These  data  (both  phase  and  ctriierence)  are  much 
better  dian  that  measured  at  the  level  groiuid  array  (refer  to  Figure  3). 

The  FWD  was  also  used  as  a  source  for  measurements  made  near  the  slope 
but  was  found  to  be  inferior  to  the  harmonic  sources.  Both  the  phase  and 
coherence  between  die  two  receivers  are  very  poor  and  inconsistent  as  shown 
in  Figure  35.  This  suggests  that  impulsive  sources  may  not  be  appropriate  to 
measure  dynamic  displacements  near  slopes. 

The  measured  displacements  were  examined  at  a  number  of  discrete  fre¬ 
quencies  between  50  and  150  Hz.  Typical  results  of  slope  measurements  are 
shown  in  Figures  36  and  37.  These  figures  present  the  real  and  imaginary 
parts  of  measured  displacements,  wrapped  and  unwrapped  phase  for  propaga¬ 
tion  down  the  slope  and  along  the  crest,  and  the  difference  between  the  phase 
in  these  two  directions. 

Comparisons  between  the  real  and  imaginary  parts  of  displacement  (refer 
to  Figures  36a  and  37a)  show  that  the  magnitudes  of  displacement  are  some¬ 
what  similar.  The  distributions  of  displacement  with  distance  are  not  very 
smooth  because  of  the  few  points  used  to  represent  the  variation.  The  ampli¬ 
tude  of  displacements  generally  increased  as  frequency  increased  but  this  is 
attributed  to  an  increase  in  force  output  of  the  vibrator  within  the  range  of 
frequencies  swept. 

An  analysis  of  the  difference  in  phase  created  by  the  slope  (refer  to  Fig¬ 
ures  36b  and  37b)  indicates  that  a  phase  lag  does  occur  for  the  sloping  (down) 
condition.  This  phase  lag  sqipears  to  be  fairly  constant  as  a  function  of  sur¬ 
face  distance,  especially  when  d|  and  dj  are  less  than  5X.  (At  much  larger 
distances,  the  difference  in  phases  has  a  constant  change  with  distance.) 
Variability  in  the  calculated  phase  difference  is  attributed  to  variability  in  the 
soil  deposit  and  sensor  distances.  These  findings  are  consistent  with  the 
results  of  numerical  iq)proximations  at  distances  greater  than  3X/4  beyond  the 
discontinuity. 

Comparisons  were  also  made  in  terms  of  the  magnitude  of  phase  velocity 
between  level  ground  and  slope  measurements  in  the  form  of  dispersion 
curves  (wavelength  versus  frequency  and  velocity  versus  depth).  Wavelength 
is  determined  from  the  slope  of  the  unwrapped  phase  relation  (versus  dis¬ 
tance).  Phase  velocity  may  then  be  calculated  using  Equation  2.7.  Wave¬ 
length  was  calculated  at  every  eighth  frequency  (increments  of  3.0  Hz) 
beginning  at  35  Hz  and  continuing  to  179  Hz  using  a  least-squares  fit  of  the 
unwTjqiped  phase  at  locations  where  the  coherence  between  the  sources  and 
receiver  was  greater  than  0.99. 

The  data  show  that  a  discontinuity  and  slope  generally  do  not  disperse 
surface  waves.  The  data,  in  terms  of  wavelength  versus  frequency,  are  shown 
.li  Figure  38  for  each  slope  site  and  direction.  Contours  of  velocity  are  also 
shown  along  with  the  best-fit  relationship  derived  from  level  ground  measure¬ 
ments.  The  wavelengths  from  the  two  locations  and  two  directions  are  very 
similar  except  when  /  >  125  Hz  where  some  scatter  is  produced  for 
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Spectrum  Cross  Power  Spectrum 
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difference  and  coherence  between  receivers  for  sloping  ground  measurements  and  impulsive  sources 
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Figure  36.  Results  for  measurements  at  Location  B  and  frequency  of  55  Hz 


Coherence  =  0.99 
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Compilation  of  results  from  slope  measurements  showing  range  in  phase  velocities  and  comparison  with  level  groutul 
measurements 


Locatioo  B  and  is  attributed  to  aliasing.  The  data  for/  <  125  Hz  generally 
fUl  within  a  band  of  phase  vdocUies  between  400  and  450  and  come 
doser  to  levd  ground  measurements  as  X  increases. 

The  data,  in  terms  of  wavdengtb  versus  phase  vdocity  (disposion  curve) 
for/  <  125  Hz,  are  shown  in  Figure  39a  along  with  the  b«t-fit  rdationship 
derived  from  levd  ground  measurements.  The  avwage  phase  change  over  the 
Imgth  of  the  array  was  calculated  using  a  least  squares  fit  algorithm.  Average 
[diase  rdations  are  shown  as  straight  lines  in  Figures  29  and  30.  These  phase 
vdocities  tend  to  frdl  within  a  narrow  band  between  400  and  460  ^s  but  are 
significantly  larger  than  the  best-fit  rdationship  for  level  ground  measurements 
and  b^in  to  deviate  more  as  X  decreases  bdow  9  ft.  Phase  vdocities  along 
the  crest  are  generally  greater  than  phase  vdocities  measured  down  the  slope. 

Phase  vdocities  measured  in  slewing  ground  were  dso  compared  with  levd 
ground  measurements  at  a  spacing  of  eight  feet.  Measured  values  of  phase 
corresponding  to  die  first  and  ninth  receiver  (di  =  10  ft  and  dj  =  18  ft)  were 
subtracted  and  the  vdocity  calculated  and  compared  with  the  best-fit  relation¬ 
ship  and  the  8-foot-spacing  data  (F4)  from  the  free  field  as  shown  in  Fig¬ 
ure  39b.  A  close  match  exists  between  phase  vdocities  measured  in  the  free 
fidd  and  down  the  slope  dthough  the  slope  vdocities  tend  to  be  slightly 
greater. 

The  phase  velocities  may  be  greater  along  the  crest  because  of  site  grading 
and  compaction  that  probably  occurred  during  devdopment  of  the  area.  The 
effort  of  field  compaction  is  likdy  to  be  greater  near  the  crest  and  less  down 
the  slope  which  would  tend  to  produce  greater  velocities  near  the  crest. 


Chapter  5  FMd  Maasuramanta 


57 


WAVELENGTH,  X  (ft)  WAVELENGTH, 
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esuits  using  least  squares  fit  of  all  data 

PHASE  VELOCITY,  C  o  (fps) 
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b.  Results  using  8-ft  spacing  and  dj/d^  =  2.3 


Rgure  39.  Dispersion  curve  from  slope  measurements  and  comparison  with  level  ground 
measurements 
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6  Summary  and  Conclusions 


The  effect  of  sloping  ground  on  the  propagation  of  surface  waves  has  been 
examined  using  a  semi-analytical  method  to  calculate  dynamic  vertical  dis¬ 
placements  produced  by  a  3-D  harmonic  load  in  planar  geosystems.  Corrobo¬ 
rative  measurements  of  vertical  dynamic  displacements  at  a  test  site  on  the 
WES  Reservation  have  also  been  reported.  The  results  of  computations  and 
field  measurements  suggest  that  surface  wave  measurements  can  be  reliably 
and  accurately  made  at  sloping  ground  sites. 

The  formulation  involves  creating  a  3-D  dynamic  stiffness  matrix  and  then 
condensing  the  components  into  an  equivalent  2-D  dynamic  stiffness  matrix. 
The  out-of-plane  loads  are  rqiresrated  by  a  Fourier  expansion  and  applied  as 
nodal  forces.  The  solution  to  the  system  of  equations  is  made  for  each  spatial 
wavenumber  and  then  the  inverse  Fourier  transform  produces  the  complex 
dynamic  displacements.  This  m^od  provides  very  accurate  values  of  static 
and  dynamic  displacements.  Brief  descriptions  of  Ae  formulation,  validation, 
and  parametric  analysis  have  been  presented. 

The  effect  of  a  slope  on  surface  wave  propagation  is  to  create  a  phase  shift 
for  waves  that  intersect  the  discontinuity.  The  phase  shift  is  generally  positive 
(phase  lag)  for  sloping  down  conditions  and  negative  (phase  lead)  for  sloping 
up  conditions.  Measurements  along  the  crest  and  base  of  the  slope  are  not 
affected  by  the  presence  of  the  slope.  In  general,  the  set-back  distance  and 
frequency  of  excitation  were  found  not  to  affect  measurements  in  any 
direction. 

Close  comparisons  were  found  to  exist  between  the  change  in  phase  down 
the  slope  and  along  the  crest  both  in  numerical  approximations  and  field  mea¬ 
surements  conduaed  at  a  site  at  WES.  Comparisons  were  also  made  between 
dispersion  curves  derived  from  numerical  approximations  and  field  measure¬ 
ments.  The  phase  velocities  near  the  slope  as  calculated  using  the  approxima¬ 
tion  method  tended  to  be  slightly  greater  than  phase  velocities  measured  in  the 
near  fidd.  Natural  variations  in  elastic  moduli  at  the  test  site  could  easily 
account  for  the  magnitude  of  difference  that  exists.  It  also  appears  that  har¬ 
monic  sources  may  be  preferable  for  measurements  in  the  vicinity  of  sloping 
ground. 

The  assumption  of  level  ground  made  for  inversion  schemes  does  not 
appear  to  be  a  very  restrictive  limitation.  The  magnitude  of  the  phase  shift 
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does  not  need  to  be  known  if  the  receivers  are  placed  sufficiently  far  away 
from  the  edge  of  die  slope.  The  results  are  independent  of  the  frequency  of 
excitation  and  set-back  distance  given  that  the  load  is  far  enough  from  the 
discontinuity. 

In  conclusion,  it  appears  that  non-destructive  surface  wave  measurements, 
including  the  current  SASW  method,  can  be  easily  adapted  to  provide  a  useful 
means  to  obtain  material  properties  in  more  complex  geosystems.  Simple 
relationships  may  be  incorporated  into  inversion  schemes  to  account  for  slope 
discontinuities. 
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